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1.  INTRODUCTION 


The  in-flight  nosetip  temperature  response  of  kinetic  energy  (KE)  projectiles  is  of  sig¬ 
nificant  interest  to  the  U.S.  Army.  The  design  of  a  new  class  of  hypersonic  KE  projectiles 
will  partially  depend  upon  a  knowledge  of  the  transient  thermal  response  due  to  free-flight 
aerodynamic  heating.  The  Aerodynamics  Branch  of  the  Weapons  Technology  Directorate. 
Army  Research  Laboratory  is  developing  and  maintaining  a  capability  to  accurately  model 
the  heating  response  of  KE  projectiles.  A  major  area  of  focus  is  the  blunt  nosetip  heating 
problem;  in  particular,  the  modeling  of  aerodynamic  surface  heat  transfer  characteristics. 

This  report  documents  recent  computational  predictions  of  nosetip  heat  transfer  rates 
for  an  M829-like  projectile  configuration.  The  M829  projectile  is  currently  launched  from  a 
120mm  gun  at  a  velocity  of  about  1.7  km/sec  (about  Mach  4.9)  at  sea-level  conditions.  The 
M829  configuration  is  serving  as  one  starting  point  for  the  design  of  hypersonic  projectile 
configurations  to  be  launched  at  still  higher  velocities.  The  heat  transfer  results  presented 
here  provide  a  model  for  boundary  conditions  in  a  separate  analysis  of  the  projectile  unsteady 
thermal  response. 

The  computational  predictions  are  made  for  free-stream  velocities  up  to  3  km/sec  (Mach 
8.8),  using  two  ditfereut  computational  approaches.  The  first  is  a  perfect  gas,  time-dependent. 
Navier-Stokes  (N-S)  finite- difference  computational  technique.  The  second  is  a  boundary- 
layer  engineering  technique,  known  as  the  ABRES  Shape  Change  Code  ( ASCC).  Comparison 
is  made  between  the  results  of  these  two  approaches  for  laminar  and  turbulent  flow  condi¬ 
tions.  As  a  basis  for  experimental  validation,  comparison  is  also  made  with  wind-tunnel  heat 
transfer  meaisurements  for  a  hemisphere-cylinder  configuration  at  Mach  6.82. 

2.  BACKGROUND  REMARKS 

Designer  concerns  about  KE  projectile  nosetip  ablation  have  recently  intensified  due  to 
the  current  interest  in  hypersonic  launch  technology.  A  major  opportunity  exists  to  apply 
and  extend  in-house  computational  fluid  dynamics  (CFD)  capabilities  in  order  to  impact  the 
design  of  the  new  class  of  hypervelocity  projectiles. 

The  complete  methodology  involves  two  major  aspects  to  be  addressed.  First,  the  in¬ 
flight  surface  heat  transfer  characteristics  must  be  determined  through  a  consideration  of  the 
projectile  flowfield.  Then,  using  a  model  which  incorporates  the  heat  transfer  characteristics, 
the  in-depth  transient  thermal  response  is  determined  by  considering  a  heat  conduction 
analysis  within  the  projectile  body.  The  long-standing  approach  is  to  consider  these  two 
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aspects  of  the  problem  separately.  This  report  specifically  addresses  the  first  aspect  only, 
focusing  on  the  heat  transfer  characteristics  of  the  blunt  projectile  nosetip. 

Among  the  important  issues  to  consider  in  the  study  of  aerodynamic  heat  transfer  is 
the  determination  of  boundary- layer  transition.  Criteria  exist  for  estimating  the  location 
and  extent  of  the  transition  region,  but  the  approach  adopted  here  is  to  compute  the  fully 
laminar  and  fully  turbulent  flow  cases  separately.  The  laminar  and  turbulent  cases  basically 
represent  the  lower  and  upper  bounds,  respectively,  of  heating  which  would  be  expected  to 
occur  in  actual  flight. 

.Although  the  Navier-Stokes  computational  technique  discussed  herein  cissumes  perfect 
gas  behavior,  it  is  emphasized  that  in-house  work  is  ongoing  to  assess  the  importance  of  real 
gas  effects  on  the  heat  transfer  characteristics  for  the  conditions  of  interest  to  the  .Army. 
The  perfect  gas  results  presented  here  will  serve  as  a  benchmark  with  which  to  compare  and 
validate  future  CFD  codes  having  perfect  gas  and  real  gas  capabilities. 

3.  COMPUTATIONAL  PREDICTIVE  APPROACHES 

Predictions  of  heat  transfer  characteristics  are  made  using  two  different  computational 
approaches.  The  first  is  a  Navier-Stokes  numerical  technique;  the  second  is  an  engineer¬ 
ing  technique  known  as  ASCC.  Comparison  is  made  between  the  results  of  these  two  ap¬ 
proaches  in  order  to  provide  a  benchmark  for  future  applications.  The  ASCC  code  is  attrac¬ 
tive  because  it  executes  quickly  on  a  typical  workstation  or  mini-computer.  Navier-Stokes 
technology,  however,  is  required  for  more  complex  problems,  such  as  those  involving  (I) 
non-axisymmetric  shapes,  (2)  flat  nosetip  shapes,  or  (3)  three-dimensional  effects  due  to 
angle-of-attack  or  geometric  irregularities. 

3.1  Unsteady  Navier-Stokes  Technique  (N-S).  Surface  heat  transfer  characteris¬ 
tics  were  computed  using  an  unsteady  (i.e.,  time-dependent)  Navier-Stokes  (N-S)  computa¬ 
tional  technique  on  a  Cray  X-MP/48  supercomputer  located  at  Aberdeen  Proving  Ground. 
This  technique,  as  reported  by  Pulliam  and  Steger  (1980)  and  Kutler,  Pulliam,  et  al  (1980), 
is  a  three-dimensional,  finite- difference,  viscous  flow  solution  procedure  for  compressible 
flow  fields.  The  N-S  technique  temporally  integrates  the  dimensionless,  transformed,  time- 
dependent,  thin-layer,  mass-averaged  Navier-Stokes  equations  in  strong  conservation  law 
form,  for  a  perfect  gas.  The  solution  vector,  Q,  is 

Q  =  {p,pu,pv,pw,e)  (1) 


2 


where  the  density  is  />;  the  Cartesian  velocity  components  are  u,  v,  and  in;  and  the  total 
energy  per  unit  volume  is  e.  The  solution  vector  is  obtained  at  each  grid  point  by  integrating 
toward  a  steady  state.  Currently,  a  set  of  inviscid,  viscous  laminar,  and  viscous  turbulent 
axisymmetric  flow  solutions  at  a  single  Mach  number  takes  on  the  order  of  an  hour  of  (^’PU 
time  to  generate  on  a  Cray  X-MP  using  a  50  x  65  grid.  Efforts  are  continuing  in-house  to 
reduce  the  time  requirements  further,  but  are  not  a  topic  for  discussion  in  this  report. 

This  Navier-Stokes  technique  has  been  applied  and  validated  in-house  (Sturek,  Guidos, 
et  al  1983;  Guidos,  Weinacht,  Dolling  1990)  for  blunt  nosetip  configurations.  A  range  of 
comparisons  has  been  made  in  the  references  between  N-S  results  and  measured  surface 
pressures,  aerodynamic  forces,  velocity  profiles,  and  skin  friction  coefficient.  Details  of  t  he 
technique,  including  the  governing  equations  and  numerical  algorithm,  can  be  traced  through 
the  aforementioned  references. 

The  local  vector  of  surface  heat  transfer  rate  per  unit  area,  q,  is  obtained  from  the 
computed  flowfield  solutions  using  Fourier’s  Law,  i.e: 

dT* 

•7  —  )  Itua//  (“) 

where  «  is  the  coefficient  of  conductivity  of  the  fluid  and  n  is  the  physical  coordinate  per¬ 
pendicular  to  the  surface.  The  vector  q  is  obtained  from  the  computed  flowfiekis  using  a 
first-order,  one-sided,  finite-difference  approximation  along  the  transformed  coordinate  which 
extends  outward  from  the  wall.  The  local  heat  transfer  rate  from  the  fluid  to  the  body  per 
unit  area,  7,  is  then  defined  as: 

q  =  -nq  (3) 

3.2  Boundary-Layer  Engineering  Technique  (ASCC).  Surface  heat  transfer  char¬ 
acteristics  were  also  computed  using  an  axisymmetric  boundary-layer  engineering  technique 
known  as  the  ABRES  Shape  Change  Code  -  ASCC  (Suchsland  1980).  The  .ASCC  code 
utilizes  a  modified  Newtonian  surface  pressure  distribution  and  inviscid  flowfield  approxima¬ 
tions  to  generate  boundary-layer  edge  flow  conditions. 

The  ASCC  code  is  of  special  interest  because  it  generates  adiabatic  wail  temperatures 
and  heat  transfer  coefficients  in  less  than  a  minute  of  CPU  time  on  a  typical  workstation  or 
mini-computer.  Such  quickness  makes  it  a  popular  tool  in  aeroballistics  design.  Comparison 
of  ASCC  results  with  Navier-Stokes  results  at  sea-level  atmospheric  flight  conditions  is  an 
important  step  in  establishing  a  benchmark  for  future  applications  of  both  codes. 
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4.  RESULTS 


4.1  Hemisphere-Cylinder  Wind-Tunnel  Case.  Hypersonic  surface  heat  transfer 
rate  measurements  for  a  hemisphere- cylinder  wind-tunnel  model  are  documented  in  a  report 
by  Crawford  and  McCauley  (1956).  The  configuration  is  a  thin- walled  hemisphere-cylinder 
of  diameter  76.84  mm,  shown  in  Figure  1.  Nominal  free-stream  conditions  are  given  as; 
Mach  number,  Moo=6.82;  Reynolds  number  based  on  body  diameter,  i?ed=1.03x  10®  and 
3.4x10®;  static  temperature,  K,  and  total  temperature,  T’o=611  K.  Real  gas  effects 

in  this  experiment  are  expected  to  be  negligible,  since  the  total  temperature  is  significantly 
below  the  threshold  which  excites  the  vibrational  modes  of  O2  or  N2.  The  flow  is  reported 
to  be  laminar,  and  the  computational  results  are  presented  here  assuming  fully  laminar  flow 
conditions. 

In  the  experiment,  the  surface  temperature  distribution  on  the  model  was  measured 
using  an  array  of  thermocouples.  Thermocouple  readings  v/ere  recorded  during  each  run 
every  3  seconds  as  the  model  temperature  increased.  The  local  heat  transfer  rate,  q,  at  each 
thermocouple  location  was  determined  through  an  analysis  which  is  described  in  the  above 
referenced  report.  In  the  report,  the  rtiCasured  heat  transfer  rates  were  then  plotted  as  a 
function  of  the  local  wall  temperature,  Ty,.  Then,  using  the  assumption  that  q  is  linear  with 
respect  to  Tu,,  linear  fits  were  determined  and  plotted.  The  differences  in  the  measured  and 
fitted  values  of  heat  transfer  rates  as  reported  suggest  a  measurement  uncertainty  on  the 
order  ±15%. 

The  Navier-Stokes  computations  were  performed  by  prescribing  a  wall  temperature  which, 
unlike  that  of  the  experiment,  is  constant  with  respect  to  both  time  and  body  location.  The 
lowest  temperature  where  the  measured  values  of  heat  transfer  rates  (in  addition  to  the  fitted 
values)  are  reported  at  all  thermocouple  locations  and  for  all  Reynolds  numbers  is  T’„,=3S5 
K  (about  63%  of  the  nominal  free-stream  total  temperature).  This  is  the  wall  temperature 
prescribed  in  the  computational  results  to  be  presented. 

Figure  2  shows  the  comparison  of  computed  and  measured  heat  transfer  rates  plotted  as 
a  function  of  axial  location  for  the  case  /?e<i=l.03x  10®.  At  the  stagnation  point  (x/rn=:0.0), 
the  computed  heat  transfer  rate  is  about  5%  higher  than  the  measured  value.  .At  x/r„=0.5 
and  x/r„=1.0  ,  the  computed  values  are  about  15%  lower  than  the  measured  values. 

Figure  3  shows  the  comparison  for  the  case  /?ed=3.4x  10®.  The  heat  transfer  rates  at 
this  Reynolds  number  are  about  60%  of  the  heat  transfer  rates  of  the  previous  figure.  The 
computed  stagnation  point  heat  transfer  rate  is  about  10%  higher  than  the  measured  value. 
At  x/r„=0.5,  the  computed  value  is  about  15%  higher  than  the  measured  value;  at  x/r„=1.0. 
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the  computed  value  is  about  10%  lower  than  the  measured  value.  The  comparisons  of  heat 
transfer  rate  shown  in  Figures  2  and  3  appear  to  be  within  the  uncertainty  of  the  measured 
data,  and  serve  to  validate  the  Navier-Stokes  computational  analysis  for  these  conditions. 

4.2  M829-Like  Configuration  Ffee-Flight  Case.  The  M829  projectile  configura¬ 
tion,  Figure  4,  is  currently  serving  as  one  starting  point  for  the  design  of  new  hypersonic 
configurations.  The  reference  diameter  is  that  of  the  cylinder,  27.05  mm.  The  standard 
nosetip  geometry,  Figure  5,  is  assumed  to  be  a  spherical  cap  of  radius  1.65  mm  tangent  to 
an  8°  cone.  Flight  conditions  are  taken  to  be  standard,  atmospheric,  sea-level  conditions. 
Results  are  shown  for  free-stream  velocities  of  1.7,  2.0,  2.5,  and  3.0  km/sec;  the  correspond¬ 
ing  Reynolds  numbers  based  on  the  reference  diameter  are  3.00  x  10®,  3.62  x  10®,  4.53  x  10**, 
and  5.40  x  10®,  respectively.  The  angle-of-attack  is  0°;  the  wall  temperature  and  free-stream 
temperature  are  both  293  K. 

Figures  6  through  9  show  the  computed  surface  heat  transfer  rates  for  velocities  of  1.7.  2.0. 
2.5,  and  3.0  km/sec,  respectively.  The  N-S  and  ASCC  solutions  are  compared  for  laminar 
and  turbulent  flow  conditions.  The  laminar  stagnation  point  heat  transfer  rates  compare 
with  a  maximum  disagreement  of  about  5%.  Downstream,  the  comparison  is  improved  to 
within  2%.  The  turbulent  stagnation  point  heat  transfer  rates  compare  to  within  about  5%. 
The  N-S  results  show  the  maximum  turbulent  heat  transfer  rates  to  be  up  to  40%  higher 
than  the  ASCC  results.  Further  downstream,  disagreement  is  also  as  high  as  40%.. 

Both  codes  show  the  maximum  turbulent  heat  transfer  rate  to  occur  downstream  of  the 
stagnation  point,  in  the  range  0.1  <  x/r„  <  0.15  (about  25°-30®  off-axis).  At  each  velocity, 
the  ASCC  results  predict  the  same  value  of  stagnation  point  heat  transfer  for  the  laminar  and 
turbulent  cases.  The  Navier-Stokes  results,  however,  predict  a  turbulent  stagnation  point 
heat  transfer  rate  about  10%  higher  than  the  laminar  rate.  In  all  cases,  the  turbulent  .ASCC 
prediction  is  bounded  by  the  turbulent  N-S  prediction  and  the  two  laminar  predictions. 

Figure  10  shows  predictions  of  laminar  heat  transfer  rate  at  2.5  km/sec  for  two  additional 
nosetip  radii  other  than  the  standard  1.65  mm  (0.065  in),  specifically  0.762  mm  (.03  in)  and 
2.54  mm  (0.1  in).  In  both  cases,  the  laminar  stagnation  point  heat  transfer  rates  between  the 
two  codes  agree  to  within  5%.  Agreement  over  the  remaining  spherical  portions  is  around 
10%,  and  agreement  on  the  cone  portions  is  within  1%.  Compared  to  the  standard  nosetip, 
the  larger  (2.54  mm  radius)  nosetip  decreases  the  laminar  stagnation  point  heat  transfer  rate 
by  about  20%.  The  smaller  (0.762  radius)  nosetip  increases  the  laminar  stagnation  point 
heat  transfer  rate  by  about  35%. 

Figure  11  shows  the  comparison  of  turbulent  heat  transfer  rate  for  the  two  additional 
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nosetip  radii.  The  relative  comparisons  between  the  two  codes  are  similar  to  those  of  the 
standard  nosetip.  The  turbulent  stagnation  point  heat  tr2msfer  rates  compare  to  within 
about  10%.  The  N-S  results  show  the  maximum  turbulent  heat  transfer  rates  to  be  about 
40%  higher  than  the  ASCC  results.  Compared  to  the  standard  nosetip,  the  larger  nosetip 
decreases  the  turbulent  stagnation  point  heat  transfer  rate  by  about  20%,  and  by  about  15% 
at  the  maximum  heating  point.  The  smaller  nosetip  increases  the  turbulent  stagnation  point 
heat  transfer  rate  by  about  40%,  and  by  about  30%  at  the  maximum  heating  point, 

5.  CONCLUSION 

Computational  predictions  have  been  presented  of  heat  transfer  characteristics  for  the 
sphere-cone  nosetip  of  an  M829-like  projectile  configuration.  Results  were  presented  for 
velocities  in  the  range  1.7  to  3.0  km/sec,  laminar  and  turbulent  flow  conditions.  Compar¬ 
isons  were  made  between  two  computational  techniques:  (1)  a  time-dependent  Navier-Stokes 
technique  and  (2)  the  ABRES  Shape  Change  Code  boundary-layer  technique. 

Comparisons  of  laminar  nosetip  heat  transfer  rate  between  the  two  codes  showed  agree¬ 
ment  within  5%  at  the  stagnation  point,  and  within  2%  further  downstream.  Comparisons  of 
turbulent  heat  transfer  rates  between  the  two  codes  showed  notable  disagreement.  The  N-S 
approach  predicted  a  10%  higher  turbulent  heat  transfer  rate  than  ASCC  at  the  stagnation 
point,  and  up  to  40%  higher  turbulent  heat  transfer  rate  at  the  point  of  maximum  heat 
transfer  rate  and  downstream.  Both  codes  showed  the  location  of  maximum  turbulent  heat 
transfer  rate  to  occur  about  25°-30'’  off-axis. 

At  each  velocity,  the  ASCC  turbulent  heat  transfer  rate  was  bounded  by  the  N-S  turbulent 
heat  transfer  rate  and  the  two  laminar  rates.  Furthermore,  results  for  two  additional  nosetip 
radii  at  2.5  km/sec  showed  comparisons  which  were  similar  to  the  standard  nosetip  radius. 

Comparison  was  also  made  between  Navier-Stokes  results  and  existing  wind-tunnel  lam¬ 
inar  heat  transfer  rate  measurements  for  a  hemisphere-cylinder  model  at  Mach  6.82.  The 
Navier-Stokes  results  at  two  different  Reynolds  numbers  compared  with  the  measured  values 
by  about  5%-10%  at  the  stagnation  point  and  within  10%-15%  elsewhere.  The  estimated 
uncertainty  of  the  measured  values  is  about  ±15%. 

The  heat  transfer  characteristics  which  have  been  presented  here  provide  a  necessary 
boundary  condition  for  subsequent  in-depth  transient  thermal  analysis  of  M829-like  projec¬ 
tile  configurations.  The  Navier-Stokes  and  ASCC  techniques  will  continue  to  be  utilized 
to  generate  heat  transfer  characteristics  for  prospective  designs.  Additionally,  the  results 
presented  here  will  serve  as  a  benchmark  for  CFD  codes  having  real  gas  capabilities. 
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Figure  1.  Illustration  of  Hemispher-Cylinder  Wind-Tunnel  Model 
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Figure  2.  Heat  Transfer  Rate,  Hemisphere-Cylinder  Wind-Tunnel  Model,  A/=6.82,  Re, 
=  1.03  X  10® 
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Figure  4.  Illustration  of  M829  Projectile 
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Figure  5.  Illustration  of  M829-Like  Nosetip 
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Figure  6.  Heat  Transfer  Rate,  M829-Like  Nosetip,  1.7  km/sec 
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Figure  7.  Heat  Transfer  Rate,  M829-Like  Nosetip,  2.0  km/sec 
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Figure  8.  Heat  Transfer  Rate,  M829-Like  Nosetip,  2.5  km/sec 
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Figure  9.  Heat  Transfer  Rate,  M829-Like  Nosetip,  3.0  km/sec 
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Figure  10.  Heat  Transfer  Rate,  2  Alternate  Nosetips,  Laminar  Flow,  2.5  km/sec 
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Figure  11.  Heat  Transfer  Rate,  2  Alternate  Nosetips,  Turbulent  Flow,  2.5  km/sec 
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LIST  OF  SYMBOLS 


d  reference  diameter 

e  total  energy  per  unit  volume 

M  Mach  number 

n  physical  coordinate  normal  to  the  wall 

Q  Vector  of  dependent  variables  in  Navier-Stokes  technique 

q  vector  of  heat  transfer  rate  per  unit  area  in  n  direction 

q  heat  transfer  rate  from  fluid  to  body  per  unit  area 

jRe<j  free-stream  Reynolds  number  based  on  reference  diameter 

r„  nosetip  radius 

T  temperature 

u,  y,  w  Cartesian  velocity  components 
V  free-stream  velocity 

X  axial  distance  from  projectile  nosetip 


Greek  Symbols 

K  coefficient  of  thermal  conductivity 

p  density 

Subscripts 

w  wall  condition 

o  total  condition 

oo  free-stream  condition 
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